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ABSTRACT 


\ 

\ 

We  present  a critical  review  of  the  fundamental  ideas 
that  underlie  the  conventional  theory  of  acoustic  detection. 
This  review  is  constructed  around  the  suggestion  that  acous- 
tic detection  apparatus  be  employed  for  collecting  informa- 
tion about  a complete  set  of  tactically  relevant  parameters 
rather  than  for  testing  hypothesis.  This  idea  is  illustrated 
by  applying  the  Woodward  & Davies  concept  of  sample-path 
information  to  the  problem  of  detection  of  signals  in  white 
noise.  Finally,  we  argue  that  within  such  an  information- 
theoretic  framework,  information  generally  disregarded  by 
conventional  procedures  can  be  extracted  from  the  received 
voltage  history.  To  obtain  this  additional  amount  of  infor- 
mation, it  is  required  that  we  possess  a far  more  detailed 
description  of  underwater  acoustics  than  what  a conventional 
theory  would  demand.  In  return,  however,  passive  range 
measurements  based  on  the  received  voltage  history  alone 
become  conceptually  feasible. 


SECTION  I 
INTRODUCTION 


THE  SEARCH  PROBLEM 

Recognizing  the  important  role  that  searching  for  submarines  plays  in  a 
major  naval  encounter,  the  Center  for  Naval  Analyses  has  initiated  a continuing 
investigation  into  our  present  capability  of  finding  submarines  by  acoustical 
means,  with  the  expectation  that  such  an  effort  will  provide  a more  accurate 
modeling  of  the  search  process  and,  perhaps,  will  suggest  to  the  Navy  better 
ways  of  performing  its  antisubmarine  warfare  task.  In  this  first  publication 
of  the  Acoustic  Detection  Project,  we  report  on  our  current  understanding  of 
the  acoustic  search  problem  and  indicate  the  direction  we  intend  pursuing  next. 

The  observation  fundamental  to  our  discussion  is  that  ocean  phenomena, 
to  the  extent  they  affect  sound-wave  propagation  through  the  ocean  environment, 
are  the  determining  factor  in  the  effectiveness  with  which  naval  forces  search 
for  an  acoustic  target.  In  fact,  far  from  being  an  unqualified  study  of  ocean 
physics,  our  investigation  is  explicitly  aimed  at  understanding  the  manner  in 
which  our  capability  to  conduct  a successful  search  is  influenced  by  the 
physical  processes  occurring  in  the  ocean.  Therefore,  a meaningful  positing 
of  the  problem  of  ocean  modeling  requires  that  we  first  embed  this  problem 
into  a larger  mathematical  framework  dealing  with  the  transformation  of 
measured  acoustic  pressure  changes  into  tactically  useful  information.  With- 
out a processing  scheme  clearly  defined,  there  exists  no  meaningful  way  of 


determining  the  type  of  description  of  ocean  phenomenology  that  is  relevant 
to  our  stated  purpose. 


r 


In  conventional  processing,  the  acoustic  pressure  changes  or  eoui valently , 
the  voltage  history  measured  at  the  output  of  the  hydrophone  system  are  consid- 
ered as  sample  values  used  to  statistically  determine  whether  or  not  a given 
signal  is  contained  within  the  received  acoustic  pressure  field.  This  puts 
conventional  processing  squarely  into  the  realm  of  statistical  hypothesis  test- 
ing and,  hence,  requires  that  the  processor  perform  a thresholding  operation 
upon  a suitably  chosen  function  of  the  voltage  history  (reference  1).  If  target 
bearing  and  frequency  spectrum  are  also  desired,  the  processor  is  allowed, 
before  the  thresholding  operation,  to  perform  the  appropriate  Fourier  transform- 
ation upon  the  received  voltage  history.  Other  parameters  of  tactical  relevance, 
notably  the  range  to  the  target,  are  distinctly  harder  to  ascertain  within  the 
conventional  framework  and,  in  fact,  modern  processors  are  not  capable  of  per- 
forming a reliable  passive  ranging  operation. 

As  we  shall  show  in  some  detail  later,  this  processing  philosophy  demands 
that  we  model  the  ocean  by  providing  a description  of  the  statistical  properties 
of  the  signal-to-noise  ratio: 


where  E represents  the  signal  energy  content  over  the  sample  length  considered 
and  N is  the  spectrum  density  of  the  noise  power. 

-2- 


We  want  to  argue  that  a more  appropriate  processing  framework  is  obtained 
if  we  identify  processing  with  the  gathering  of  information  about  tactically 
relevant  parameters  rather  than  with  statistical  hypothesis  testing.  Clearly, 
to  render  this  idea  mathematically  meaningful,  a precise  definition  of  what  is 
mean  by  information  must  be  provided  and  we  shall  do  so  presently.  It  will  be- 
come apparent  from  the  definition  given  to  the  concept  of  information,  that  the 
amount  of  information  one  has  about  a given  parameter  is  a monotonic  function 
of  the  statistic  a conventional  processor  would  evaluate.  Hence,  in  conventional 
processing,  information  that  comes  below  the  chosen  threshold  is  entirely  dis- 
regarded, conventional  gear  appears  to  be  operating  with  less  information  than 
is  actually  available  in  the  random  voltage  history.  Furthermore , since  withii 
this  framework  one  would  evaluate  information  about  all  tactical  parameters  alike, 
passive  ranging  is  no  longer  conceptually  different  from  any  other  part  of  the 
detection  process. 

If  we  allow  this  information-theoretic  framework  of  signal  processing  to 
determine  what  features  of  the  ocean  environment  should  be  modeled,  we  find  that 
we  now  need  the  probability  density  function  of  the  sample  voltage  measurement 
conditioned  upon  the  value  taken  by  the  tactical  parameter  of  interest;  this  is 
a fare  more  detailed  description  of  the  ocean  than  that  required  by  the  conven- 
tional philosophy,  which  is  as  it  should  be  for  we  expect  this  processing  pro- 
cedure to  provide  us  with  more  information  about  relevant  parameters  than  the 
conventional  processors  do. 


Due  to  the  complexity  of  the  problem  at  hand,  we  only  claim  to  have  argued 
that  the  information-theoretic  approach  is  reasonable.  We  hope  to  show  later 
that  it  is  also  useful  in  that  it  provides  the  searchers  with  an  increased  capa- 
bility to  find  the  target  submarine. 

THE  CONCEPT  OF  INFORMATION 

In  what  follows,  we  shall  take  the  position  that  the  only  element  of  com- 
munication between  the  searcher  and  the  submarine  he  searches  for  is  the  random 
voltage  history  that  obtains  at  the  outlet  of  his  hydrophone  system.  Depending 
on  how  he  chooses  to  look  at  it,  this  voltage  history  will  provide  him  with 
information  concerning  various  parameters  of  tactical  usefulness  to  his  search 
effort.  But  before  we  discuss  the  details  of  how  this  is  done,  it  behooves  us 
to  specify,  and  subsequently  quantify,  the  concept  of  information  we  intend  to 
employ  here.  Following  the  1952  work  of  Woodward  & Davies  (reference  2),  let 
5 represent  the  tactical  parameter  of  interest,  Vg^  = {v(x),  re[0,t]}  the 
voltage  history  that  the  searcher  has  measured  over  the  [0,t]  time  interval, 
and  U(?jVgt)  the  searcher's  state  of  knowledge  at  time  t concerning  a par- 
ticular value  l given  the  specific  voltage  sample  path  Vg^  . The  conventional, 
and  perhaps  the  most  convenient  representation  for  one's  state  of  knowledge  is 
obtained  if  one  is  willing  to  consider  5 a random  variable  and  U(£|Vg*)  the 
corresponding  conditional  probability  density  given  the  measurement  Vgt  . 

When  a measurement  has  been  performed,  values  of  the  tactical  parameter 
that  would  imply  a high  probability  for  havinq  measured  Vq*  will  appear  to  the 
observer  more  believable  than  other  possible  values.  Due  to  distortion  by  the 


-4- 


random  medium  through  which  the  message  has  propagated,  this  reevaluation  of 
beliefs  will  not  lead  to  an  a posteriori  state  of  knowledge  entirely  concentrated 
on  a single  value  of  £ . Rather,  the  searcher's  state  of  knowledge  concerning 
the  value  of  some  tactical  parameter  after  having  read  the  voltage  history 
will  merely  differ  from  his  a priori  state  of  knowledge  U(f|V^)  . This, 

according  to  Woodward  & Davies,  is  the  essence  of  the  acoustic  detection  process. 


It  seems  intuitively  reasonable  that  the  amount  of  information  contained 
within  a measurement  Vq1  should  be  related  to  the  change  that  measurement  pro- 
duced in  the  observer's  state  of  knowledge.  To  quantify  this  notion.  Woodward 
& Davies  require  the  following  statements  to  be  true: 


• If  2 successive  voltage  readings  are  employed  to  learn  about 
the  value  of  some  tactical  parameter,  and  the  qbserver  regards 
his  a posteriori  state  of  knowledge  after  the  first  reading  as 
the  a priori  one  before  the  second,  the  total  gain  in  informa- 
tion concerning  the  value  of  the  parameter  is  equal  to  the  sum 
of  the  gains  from  each  reading. 

o If  a voltage  reading  is  received  and  used  to  learn  about  the 
value  of  2 independent  parameters,  the  total  gain  in  informa- 
tion concerning  them  is  the  sum  of  the  gains  when  each  param- 
eter is  considered  separately. 


From  these  2 axioms,  Woodward  & Davies  develop  the  whole  mathematical 
theory.  Their  argument  proceeds  as  follows.  If  upon  receiving  a reading  V ^ 
the  observer  is  certain  that  the  value  of  the  tactical  parameter  is  £ , the 

amount  of  information  contained  in  that  reading  depends  only  on  the  observer's 
initial  state  of  knowledge  IJ(£|Vq^)  . The  first  axiom  then  implies  that  the 
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amount  of  information  contained  in  a voltage  reading  which  leaves  the  observer 
with  some  uncertainty  as  to  the  actual  value  of  £ is  a function  of  the  a 
priori  and  a posteriori  states  of  knowledge  about  that  value  alone  and  may  be 
written  in  the  form  0[U(||Vq°)  , U(£|V(-^)]  . With  this  notation,  the  first 

axiom  states: 

o[u(i|v0°),  u(i|v0t 1 )]  + J[u(i|v0tl),  u(5|vtlt2)] 

= J[u(i|v0°),  uUiv0t2)]  . 

Woodward  & Davies  show  that  to  satisfy  this  identify,  J must  be  of  the  general 
form 

0(a,£5)  = j (a)  - j ( 6)  . 

In  order  to  determine  the  functional  form  of  j , the  second  axiom  is  used. 

Here  the  2 independent  parameters  may  be  denoted  by  % and  o , and  the 
corresponding  received  voltage  indication  again  by  . If  the  state  of 

knowledge  concerning  the  joint  entity  of  2 independent  tactical  parameters  is 
taken  to  be  the  product  of  the  states  of  knowledge  concerning  them  separately, 
the  second  axiom  gives 


j[U(i|V0°)  U( n | VR° ) ] - j[U(c| Vq1)  U(n| V0t)] 

S j[(i!V0°)]  - JtU(i|V0t)3  + j[U(n|V00)]  - j[u(njv0t)]  . 

From  this  identity,  j(U)  must  be  of  the  form 
j(U)  = A log  U + B 

where  A and  B are  constants  which  may  be  chosen  arbitrarily.  Thus, 

„ ti  . to  U( II v t2) 

J[U(5|V0  '),  U(c|VQ  -)]  = log  ~rr-  • 

U(c|V0  ') 

where  A was  so  chosen  as  to  make  an  increase  in  the  observer's  state  of 

knowledge  represent  positive  information.  This  result  of  Woodward  & Davies 

is  the  expression  we  shall  use  in  this  paper  to  represent  the  quantity  of 

t2 

information  gained  from  Vfj  concerning  the  value  of  the  tactical  parameter 

•* . 1 2 

C . As  might  be  expected,  it  is  equal  to  zero  if,  and  only  if  UU|VQ  ) 

ti  t2 

and  U(s|Vq  ) are  identical,  i.e.,  when  the  communication  Vtl  leaves  the 

observer's  state  of  knowledge  about  l completely  unchanged. 

SUMMARY 

In  section  II,  we  shall  briefly  discuss  those  factors  in  the  submarine 
search  problem  that  are  responsible  for  the  random  nature  of  the  voltage 
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received,  and  shall,  hence,  mention  the  major  aspects  of  underwater  acoustic 
propagation.  In  section  III  we  present  a number  of  various  ways  in  which  the 
received  voltage  can  be  manipulated  and  we  evaluate  the  amount  of  information 
that  can  thus  be  obtained.  In  section  IV,  we  give  a general  description  of 
how  conventional  modeling  employs  this  information  to  evaluate  Navy's  tactical 
performance.  We  conclude  this  paper  with  a critical  evaluation  of  conventional 
models  and  delineate  some  research  topics  we  intend  addressing  within  the 
Acoustic  Detection  Project. 


SECTION  II 

THL  RANDOMNESS  OF  THE  OCEAN  ENVIRONMENT 


As  already  outlined  in  the  preceding  section,  the  main  input  to  the 
acoustic  detection  process  is  the  random  voltage  measured  at  the  output  of 
the  observer's  hydrophone  system.  We  have  also  argued  that  the  essence  of 
acoustic  detection  is  the  change  that  this  random  voltage  continuously  induces 
in  the  observer's  state  of  knowledge  regarding  any  tactically  relevant  param- 
eter. That  this  state  of  knowledge  should  be  imperfect  is  a consequence  of 
the  random  nature  of  the  received  voltage  and  the  finite  nature  of  the  obser- 
vation interval.  It  is  therefore  essential  that  we  attempt  to  understand  the 
physical  reasons  for  this  randomness. 

Possibly  the  most  important  contributor  to  the  random  nature  of  the 
acoustic  voltage  is  the  ocean  medium  through  which  the  submarine  signal  must 
propagate  to  arrive  ^ the  observer's  hydrophone  system.  Indeed,  the  physical 
properties  .if  the  ocean  at  every  point  within  it  undergo  irregular  fluctuations. 
Similarly,  fhe  ph/sical  properties  of  the  ocean  at  different  spatial  points 
within  it,  but  at  the  same  instant  in  time  differ  from  one  another  in  a random 
fashion.  In  particular,  since  the  acoustic  index  of  refraction  of  the  ocean 
is  a function  of  uch  physical  properties  as  temperature  and  salinity,  one  can 
take  the  view  point  that  the  refractive  index  is  a random  field.  There  are 
many  different  mechanisms  that  have  been  proposed  for  modeling  the  random  prop- 
erties of  the  acoustic  refraction  index  of  the  ocean  (references  3,  4,  and  5). 


It  appears  unlikely,  however,  that  any  one  such  model  will  eventually  account 
for  all  of  the  observed  features  of  acoustic  fluctuations  in  the  ocean.  Rather, 
we  must  rely  in  practice  upon  a mixture  of  limited  models,  in  which  one  mechanism 
might  dominate  over  the  others  under  specified  conditions,  but  none  suffices  to 
account  for  underwater  acoustic  behavior  under  all  conditions. 

The  conventionally  accepted  procedure  for  translating  the  statistical 
properties  of  the  random  ocean  into  the  implied  statistical  properties  of  an 
acoustic  pressure  field  supported  by  such  a medium  is  to  use  the  v/ave  equa- 
tion with  suitably  designed  boundary  conditions.  Quite  often,  these  boundary 
conditions  themselves  are  random  which  only  aggravates  the  already  formidable 
task  of  finding  solutions  to  the  wave  equation.  Due  to  the  mathematical  com- 
plexity of  this  approach,  a number  of  approximate  solutions  to  the  wave  equa- 
tion under  simplifying  assumptions  about  the  acoustic  process  have  been 
developed  and  used  over  the  years  (references  6,  7,  and  8). 

It  is  not  our  intention  to  review,  at  this  time,  progress  in  the  field 
of  acoustic  propagation  in  random  media  in  any  detail;  the  complexity  of  the 
problem,  as  well  as  the  amount  of  professional  literature  that  has  been  written 
to  date,  precludes  such  an  attempt.  We  only  wish  to  observe  that  whatever  the 
accepted  theory  of  underwater  sound  propagation,  it  should  provide  the  proba- 
bility distribution  of  the  acoustic  voltage  history  measured  at  the  hydrophone 
output  in  terms  of  all  those  parameters  that  are  deemed  necessary  to  specify 
such  distribution  entirely.  Since  a parameter  that  does  not  in  any  way  affect 


- 
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this  probability  cannot  be  estimated  by  voltage  measurements,  the  set  of  param- 
eters necessary  to  specify  the  voltage  probability  distribution  will  be  referred 
to  as  a complete  set  of  tactically  relevant  parameters.  We  shall  assume  in  what 
follows  that  such  a probability  distribution  is  somehow  made  available  to  us, 
and  proceed  to  evaluate  the  amount  of  information  regarding  tactically  relevant 
parameters  that  is  contained  within  a sample  voltage  measurement. 
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SECTION  III 

[ THE  INFORMATION-THEORETIC  APPROACH  TO  DETECTION 

[ THE  EXISTENCE  PARAMETER 

The  most  straightforward,  although  not  necessarily  the  most  complete, 
application  of  Woodward  & Davies'  concept  of  sarnple-path  information  to  prob- 
lems of  signal  detection  appears  in  a 1952  paper  by  Davies  (reference  9).  In 
that  work,  the  important  tactical  parameter  is  taken  to  be  a 2-valued  existence 
parameter,  e , designed  to  represent  either  the  presence,  e = YES  , or  non- 
presence, e = NO  , of  a specified  acoustic  signal  in  the  received  voltage 
I history.  The  observer’s  state  of  knowledge  concerning  this  parameter  is  repre- 

sented by  the  conditional  probability  p ( e j Vq c ) that  the  signal  be  present 
within  the  received  voltage  history  given  the  sample  path  . To  evaluate 

the  time  dependence  of  this  probability  along  a voltage  sample  path,  the  theory 
relates  it  via  Bayes'  rule  to  the  conditional  probability  that  the  sample  path 
under  consideration  would  be  obtained  given  a particular  value  of  the  existence 
parameter, 


pUIV)  = — 


p(e|V0°)  P(V0t|a) 


p(YESjV0U)  p(VQt( YES)  + p(N0|V0°)  p(V()t|N0) 


(1) 


where  th$  probability  p(e|Vq  ) represents  the  observer's  a priori  state  of 
knowledge  regarding  the  value  of  the  existence  parameter  e • Therefore,  the 
amount  of  information  regarding  the  presence  of  a target  signal  within  the 


-12- 


J 


is  given,  in  accordance  with  the  definition  provided  in 


voltage  history  \l^ 
the  previous  section,  by 


o t pM^) 

J c p ( e 1 VQ  ),  pU|V0  )]  = log  jp  > 

P(e|v0  ) 


and  using  equation  (1), 


J[p(e|V0°),p(£|V0t)]  = log 


p(y  10 

p(YES| V0°)p(V0t[ YES)  + p(NO| V0°)p(V0t|NO) 


Equation  (1)  mathematically  specifies  the  ideal  receiver  for  estimating 
the  presence  of  a signal  within  the  wave  form  . Within  an  information 

theoretic  approach  to  detection,  all  chat  can  be  reasonably  demanded  of  a 
receiver  is  that  it  shall  enable  an  observer  to  determine  probabilities  that 
each  possible  value  of  the  parameter  of  interest  is  the  true  one.  If  the 
receiver  actually  computes  these  probabil ities,  no  further  problem  of  inter- 
preting the  content  of  the  wave  form  regarding  the  desired  tactical  parameter 
remains. 


THE  CASE  OF  A DETERMINISTIC  SIGNAL 

In  order  to  apply  the  foregoing  theory  to  a practical  problem,  the  illus- 
trative value  of  which  might  be  quite  welcome  at  this  point,  it  is  necessary  to 
consider  how  the  a posteriori  state  of  knowledge  p(£|VQt)  is  constructed 
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from  the  received  communication  L . We  shall  presume  that  the  a priori 
state  of  knowledge  p(tjV()^)  is  given  and  the  a posteriori  distribution  may 
then  be  obtained  directly  from  equation  (1)  once  p ( t e ) ^as  ^een  evaluated. 

The  conditional  probability  p ( | YES ) represents  the  unpredictable 
nature  of  the  received  wave  form  when  the  target  signal  is  known  to  be 

present  and,  hence,  describes  the  effect  various  random  disturbances  present 
in  the  reception  process  have  upon  the  shape  of  the  signal  wave  form.  In 
general,  these  random  disturbances  can  be  classified  either  as  distortions  in 
the  signal  wave  form  that  are  mainly  due  to  randomness  in  the  propagation 
medium,  or  as  interferences  from  sound  sources  other  than  the  target  such  as 
ambient  sound  producers,  thermal  noise  in  the  receiving  electronic  apparatus 
and  altogether  everything  else  that  might  contribute  to  the  conditional  proba- 
bility p(V(jt|N0)  . 

In  the  following  example,  signal  wave  form  distortions  will  not  be  con- 
sidered. Rather,  we  shall  assume  that  the  signal  wave  form  Snt  is  exactly 
known  and  that  the  voltage  history  received  in  the  presence  of  a target  differs 
from  S0l  only  in  that  a random  wave  form,  generally  referred  to  as  noise,  has 
been  added  to  it.  The  large  variety  of  contributory  phenomena  make  the  noise 
wave  form  and  its  statistical  properties  very  difficult  to  describe  mathematically. 
To  make  our  example  manageable,  we  shall  represent  the  noise  by  a stationary, 
band-limited,  white,  Gaussian  random  process.  Under  these  conditions,  the  evalua- 
tion of  p(V(.t|c)  is  fairly  straightforward  and  will  now  be  described. 


A simple  way  of  formulating  the  statistical  properties  of  the  noise  process 
is  by  means  of  wave  form  sampling  process.  Sampling  analysis  rests  on  a well- 
known  mathematical  theorem  (reference  10)  that  if  a function  of  time  F(t)  con- 
tains no  frequencies  greater  than  W , then  over  any  finite  observation  time  1 


2WT-1 
F ( t)  = £ 

j=0 


7 j\ 

sin  it  (2WT 

- j) 

\2W  ) 

2WT  sin  - 

dr) 

The  importance  of  this  identity  is  that  it  enables  a continuous  function 
of  time  to  be  specified  uniquely  in  terms  of  sample  values  at  intervals  ’ 

where  W is  an  arbitrary  frequency  greater  than  any  which  occurs  in  the  frequency 
spectrum  of  F(t)  . 

The  normality  of  the  assumed  noise  process  implies  that  at  the  output  of 
the  filter,  each  sampling  value,  n(jvjj  ’ or  more  briefly  n j , has  by  definition 
the  Gaussian  probability  density 


where  o2  is  the  mean  square  value  of  n(t)  . It  is  not  difficult  to  show 
that  these  noise  samples  are  statistically  independent.  Indeed,  under  the 
stationarity  assumption,  the  correlation  function  of  the  noise  process 


4>(t,t+-r)  = (n(t)  n(t+-r)^ 
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where  brackets  indicate  that  the  average  over  an  ensemble  is  to  be  taken  of  the 


quantity  within,  becomes  a function  of  -t  alone  (reference  11).  Its  Fourier 
transform  with  respect  to  x , 


4>(t) 


r dv  e~i2lIVT  *(v) 


(3) 


represents  the  noise  power  spectrum  (reference  11)  and,  hence,  with  our  whiteness 
assumption,  can  be  written  as 


(1»0  M<« 

$(v)  = < 

( 0 otherwise 


Performing  the  integration  indicated  in  equation  (3), 


sin  2-rrWx 


<p(t)  = NqW  2irWx 


and  noise  samples  taken  ^ apart  are  indeed  statistically  uncorrelated.  There- 
fore, the  joint  probability  density  of  a set  of  samples  is  the  product  of  each 
separate  distribution.  Since  the  samples  determine  the  wave  form,  this  product 
also  gives  the  probability  density  for  the  wave  form  itself.  With 


.2  s 


*(0)  = NqW  = N , 


we  thus  have. 


2Wt 


p ( v 1 NO ) 


i-ur  £ VJ2  • 


2 N / 


J=1 
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where  v represents  che  set  of  samples  measured  and  corresponds  to  a noise  read- 
ing of  length  t . Using  the  orthonormal  it}'  of  the  base  functions  employed  in 
the  sampling  theorem,  we  can  show  that 

2Wt  t 

v.2  = 2W  | ot  v2(t)  , 

j=l  J O' 

which  allows  the  wave  form  probability  density  to  be  written  in  a form  that  no 
longer  makes  any  explicit  reference  to  the  sampling  procedure: 

.t 

-1/Nn  I dx  v2(t) 
p(V0  1 NO)  = ke  O' 

Using  similar  techniques,  we  can  show  that  if  the  signal  is  also  present 
in  the  received  voltage  history,  the  probability  distribution  of  the  resultant 
wave  form  is: 

ft  >2 

-1/No  | dt(v(x)-s(x) 1 

p(V0t|YES)  --  ke  (T 

and,  therefore,  the  amount  of  information  contained  within  the  voltage  history 
regarding  the  presence  of  a signal  wave  form  becomes: 

d[p(YES|V0°),  p(YES|V0t)] 


-l/N0oft  dx(v(x)-s(x)) 


To  simplify  this  expression,  a dimensionless  energy  parameter  P(t)  , is 
defined  by  the  equation: 


“(t)  = 


2E 


N0  ’ 


where: 


E = I dr  s2( t)  , 
0 


is  the  signal  energy  received  during  the  interval  of  observation  t . Equation 
(4)  may  then  be  written 


J[p(VES|V0°),  p(VES|V0t)] 


2/N0  fl  dr  v(t)  s(t) 


= log- 


p(VE5|V0°)  e V 


2/Nn  dt  v ( t)  s(t)  0 p2(t)/2 

U „ J _L  r,  / M n v u \ o ' L ' ' c 


+ P(N0 | VQ  ) e 


(5) 


Notice  that  the  amount  of  information  regarding  the  presence  of  a target 
signal  depends  upon  the  measured  voltage  history  only  through  the  correlation 
integral , 


q(t)  = I dt  v ( t ) s(t) 
O' 
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Since  this  quantity  will  tend  to  be  greatest  when  the  received  wave  form  actually 
contains  the  given  target  signal,  the  information  regarding  the  signal's  presence 
within  the  received  wave  form  will  also  be  greatest  when  the  received  voltage 
history  contains  the  signal. 


The  form  that  the  receiver  should  take  in  the  case  just  considered  follows 
from  the  preceeding  theory.  Specifically,  equation  (5)  shows  that  in  order  to 
obtain  the  relevant  information,  the  received  wave  form  should  be  passed  through 
a filter  whose  response  to  a unit  imp! use  is 


g(t)  = 


j s(t-x) 

I 0 


0<T<t 

otherwise 


where,  as  before,  t is  the  period  of  observation.  Indeed,  if  the  signal  is 
completely  known,  the  output  of  the  filter  at  the  end  of  the  period  of  observa- 
tion yields  a voltage  proportional  to  the  correlation  integral  q(t)  . There- 
upon, the  amount  of  information  the  observer  has  gathered  follows  from  the  simple 
non-singular  transformation  of  the  filter  output  indicated  by  equation  (5). 
Finally,  from  equations  (1)  and  (2),  the  observer's  state  of  knowledge  regarding 
the  presence  of  the  signal  within  the  received  voltage  history  is  given  by: 

p( YES t VQt)  = p( YES | V0°)  eJ[P(YESiv0°)»  P< YES i VO1)  1 . 

Needless  to  say,  we  can  develop  the  theory  along  similar  lines  to  encompass 
the  calculation  of  the  amount  of  information  regarding  any  other  parameter  of 
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tactical  value,  and  thus  ascertain  the  form  of  the  corresponding  optimal  receiver. 
Also,  we  may  attempt  to  apply  the  approach  we  have  taken  above  to  a case  with  a 
somewhat  more  realistic  description  of  the  noise  process,  and  it  is  to  this 
latter  problem  that  we  turn  our  attention  first. 

THE  CASE  OF  A STOCHASTIC  SIGNAL 

If  applications  of  the  framework  described  above  to  the  problem  of  radar 
reception  were  our  purpose,  the  correlation  receiver  would  not  be  such  a bad 
approximation.  In  sonar  detection,  however,  the  fundamental  assumptions  made 
in  studying  the  deterministic  signal  case  are  not  acceptable.  The  most  impor- 
tant difficulty  is  due  to  the  fact  that  the  signal,  when  present,  has  a stochastic 
character,  and  the  neglect  of  signal  distortion  due  to  the  propagation  medium  is 
no  longer  a tenable  position. 

To  illustrate  the  manner  in  which  we  might  extend  the  theory  to  accommodate 
stochastic  signals,  we  adopt  a particular  model  (reference  12)  of  the  signals 
and  of  the  noise  in  which  they  are  to  be  received,  and  attempt  to  evaluate  the 
information  that  the  received  voltage  history  contains  concerning  the  existence 
parameter.  The  signal  to  be  detected  is  taken  as  a segment  of  a narrowband  wave 
with  a complex,  random  envelope  Z(t)  and  a carrier  frequency  fi  , 

s(-t)  = Re  Z(t)  e^T  . 


Here,  only  the  real  part  of  the  complex  wave  is  to  be  taken  as  the  signal  wave 
form.  Of  the  random  envelope,  we  only  know  that  its  components  are  Gaussian 
random  processes  with  given  complex  correlation  functions 


t ( T-J  ,t2)  - J (z(x-])  Z (t2)) 


where  the  asterisk  notation  represents  complex  conjugation  of  the  quantity  it 
modifies.  The  input  to  the  receiver  system  always  contains  Gaussian  noise  to 
which  the  signal,  when  it  occurs,  is  added.  We  shall  assume  for  the  sake  of 
computational  simplicity,  that  the  background  noise  is  white,  providing  for  the 
correlation  function  of  the  noise  process 


sin  27rW(T2-T-|) 
2ttW(t2-t-]  ) 


Since  the  actual  value  of  the  noise  bandwidth  W does  not  appear  in  the 
result,  we  take  the  limit  as  this  bandwidth  approaches  infinity  and  thus  have 
a Dirac  delta  correlation  function 


n(T1)  n (t2))  = Nq«(t2-t1) 


Therefore,  since  the  signal  and  the  noise  are  statistically  independent,  the 
complex  correlation  function  of  the  received  voltage  envelope  when  both  are 
present  is 

v (t2))  = ^ ’ T2 ) + N06(t2-t-|)  . 


To  proceed,  we  might  use,  as  we  have  before,  the  values  v(j/2W)  of  the 
voltage  history  at  propitiously  chosen  time  points.  However,  these  values  are 
no  longer  statistically  independent  of  each  other;  and  because  of  the  form  of 
their  joint  probability  density  function,  the  passage  to  the  continuum  is  some- 
what more  difficult  to  analyze  than  it  was  in  the  example  described  before. 

It  will  be  simpler,  rather,  to  find  a set  of  quantities  Vj  that  are  uncorrelated 
but  can  be  generated  by  linear  operations  on  the  voltage  wave  form.  These 
quantities  will  be  the  coefficients  of  an  expansion  of  the  received  voltage 
v ( t ) in  a particular  kind  of  Fourier  series  (reference  13). 

We  accomplish  such  by  writing  for  the  received  voltage, 

v(t)  = £2  V,f,(-r)  > 

j 

where  the  set  of  functions  f-(T)  is  orthonormal  with  respect  to  the  time 

interval  0<t<t  , 

„/'  fj*{T)  fk<’>  ‘ • 

Then,  it  is  not  difficult  to  show  that  the  expansion  coefficients  in  equation 
(7)  can  be  computed  from 
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The  expansion  coefficients  are  therefore  linearly  related  to  the  input  voltage 
v ( t ) , and  the  jth  one  could  be  generated  by  passing  the  voltage  through  a 

linear  filter  matched  to  a signal  of  the  form  fj(0 


Since  the  voltage  history  is  random,  so  are  the  expansion  coefficients  v^ 
Their  covariances  are  given  by  the  following  equation. 


(vjvk*)  = J dTl  dT2  fj*(Tl}  W (v(t1)  v*(t2))  * 

' 1 1 


and  using  equation  (6), 


Vk  ) 5 N0  fjk 


ft  * 

jk  + 0J  dTl  dT2  ♦(t1*t2)  fj  (t1}  Vt2} 


(8) 


Our  goal  will  be  achieved  if  we  can  find  a set  of  orthonormal  functions  such 
that  this  expression  vanishes  whenever  j and  k are  different.  This  will 

be  so,  if  for  each  function  f-(x)  separately  we  require 

J 


j dx  fj(-r)  - Aj  fj('c')  , (0<x"<t) 


(9) 


where  the  constants  ,\j  have  yet  to  be  determined.  Substituting  equation  (9) 
into  the  covariance  equation  (8),  we  find 


(vj  \*>  * <Y"o)sjk 
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so  that  the  expansion  coefficients  are  indeed  uncorrelated,  and  their  variances 
are  simply  related  to  the  eigen  values  (reference  14)  A of  the  integral 
equation  (9). 


A . 


+ No 


Hence,  the  joint  probability  density  of  the  entire  set  of  coefficients  is 
again  the  product  of  each  separate  distribution  and,  since  the  coefficients 
determine  the  wave  form,  this  product  gives  the  probability  density  for  the  wave 
form  itself.  But,  because  the  voltage  coefficients  v.  are  formed  by  linear 
operations  on  the  Gaussian  process  v(t)  , they  must  be  Gaussian  random  variables, 
and  we  finally  obtain, 


M 

P(v 0tiYEs)  = iim  n 

M -x  “ j = 1 


M 

-Zkl 

J=t  J 


2/2(Aj+NQ) 


Furthermore,  if  noise  alone  were  present,  the  probability  density  for  the  voltage 
history  would  trivially  be  given  by 

1 M 2 

, / , \"  -^o  ^|vjl 

p(»  |N°)-n™/— -)  e J-l 

and,  therefore,  the  amount  of  information  concerning  the  presence  of  the  target 
signal  within  the  voltage  history  already  viewed  by  the  observer  becomes 
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J[p(VES !V0°) , p(VES|V0t)] 


= log 


V2n0  jV,  \jlv//(yV 


on,  , * ]/2%  | /,1'jl  /(yN0) 

e<VESlvo>  J?i  J 


, (10) 


+ p(NO} V0°) 


provided 

.t  2 

^ j dx^  dtjUU-j  ,T^)  | <® 

to  ensure  that  the  limit  as  the  number  of  coefficients  t;ends  to  infinity  does 
indeed  exist  (reference  15). 


Before  we  can  attempt  to  determine  the  design  of  the  optimum  detector  sug- 
gested by  equation  (10),  we  must  cast  this  equation  into  a more  suitable  form. 
Thus,  employing  the  equation  defining  the  voltage  coefficient  Vj  , 


dT-|  dT^  v ( t*|  ) h(t^  ,T2>1/N)  vf-r^) 


where  we  have  defined 


h<vvu>  - £ fjO,)  y<uy) 


on 


Furthermore,  the  infinite  product  in  equation  (10), 


D(l/N0)  e n (l+^j/N0)  • 
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can  also  be  evaluated  in  terms  of  the  function  defined  in  equation  (11).  Indeed, 
taking  the  logarithm  of  D(1/NQ)  , we  have 


lnD(l/Nn)  = £ ln(l+A./Nn)  = L J 
0 j-1  J 0 j=l  0 


VNr 


du 


+ A . u 


(12) 


But,  setting  f|  = T2  in  equation  (11)  and  integrating  from  0 to  t , we 
obtain,  upon  using  the  normalization  condition  for  the  f.  functions,  the 

J 

integral  on  the  right-hand  side  of  equation  (12),  Therefore  (reference  16), 


fl/N  0 fL 

J du  / 


dx  h(x,x;u) 


D(l/N0)  = 


and 

J[p(YES|V0°),  p(YES|V0t)] 

1/2Nq2  f dx-|  dx2  v ( x-j ) h(x1  ,x2;1/Nq)  v(x?) 

_ -]0g ^ — — _ (13) 

1/2N()  [ dx-j  dx2  v*(x^  h(x^  ,x2,1/Nq)  v(x2) 

p(YESlV0°)e  0 ' + P(N0|V0°)  D(l/N0) 

indicating  that  a knowledge  of  h(x-|,x2;u)  is  sufficient  to  provide  the  informa- 
tion contained  within  the  given  sample  path  Vg  . If  we  multiply  both  sides  of 
equation  (11)  by  u<t>(x.x-,)  and  integrate  over  the  range  0<x1<x  , we  get,  using 
the  integral  equation  (9) 
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r 


I dr-  :'t  / hi-.,,x,;u)  - 'L.  f - ( x)  f.  < t)  X.  /(l+>,./u) 

1 1 - i=1  J J 4 u J 


O' 


which,  upon  addition  t nation  (11),  provides 


h(iT-t0;u'  u > dx  4>( t,  , t ) h(x,T?;u)  = Z..  A , ^(ii)  ^ (to) 
1 ' ,V  1 ^ v-l  J ■’  ' 


Finally,  since 


j- ‘ ° 


as  can  be  easily  shown,  we  find  that  h(x.|,xp;u)  is  the  solution  of  the 
following  in*egral  egu.  non: 


h(T^,^2»U)  + U / dx  <J>(  T-J  » t)  h(T  » Xp  s U ) = <f  ( X -j  , Tp ) * 0<(Ti»Tp)<t 


(14 


We  are  now  ready  to  specify  the  optimum  receiver  implied  by  equation  (13; 
. r.  tors  of  the  solution  to  the  integral  equation  (14)  To  do  sc,  let 


W(Tl)  - , 1 in,  h(Tl  ,x0;l/iJ)  v(T?)  , 

I rj  - • ^ 


and  tin 


iS  write  the  < 'rents  in  equation  (13)  as 
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J 


,.t 


1/2N 


lll 


★ 

v (ti 


) h(- 


,t2;1/n0) 


v(t2)  - 


1/N 


!l 


dr. 


v (t. 


) W(t1 ) 


For  each  value  of  , W(t-j)  is  a weighted  average  of  the  input  voltage  history 
v(t)  that  has  been  viewed  previously  to  the  time  t-j  . It  can  be  found  by 
passing  the  voltage  received  through  a time- variable  linear  filter  having  a 
narrowband  impulse  response 


9 


(t) 

1 


j h(  x | ,t-|-t,1/N) 

I 0 


0<T<T-| 

otherwise 


1 

The  output  of  this  filter  at  time  x is  Re  W(T)e  . Therefore,  if  we 
now  multiply  this  at  each  instant  by  the  input  Re  v(t)einT  and  thereupon  integrate 
the  result  with  a low-frequency  filter,  the  output  at  time  t is  essentially 
equal  to  the  exponent  in  equation  (13),  and  this  system  evaluates  the  amount  of 
information  contained  within  the  received  voltage  history  regarding  the  presence 
of  a Gaussian  narrowband  signal. 


PASSIVE  RANGING 


We  have  applied  the  sample-path  information  concept 
to  the  estimation  of  the  existence  parameter  under  varied 
underlying  statistical  properties  of  the  voltage  process, 
no  objection  to  ev tending  the  application  of  this  concept 
culatioo  of  thu  amount  of  information  regarding  any  other 
value.  In  particular,  one  wants  to  estimate  the  position 


of  Woodward  & Davies 
assumptions  about  the 
There  exists,  however, 
to  encompass  the  cal- 
parameter  of  tactical 
of  the  target  with 


r 


respect  to  some  specified  system  of  reference,  by  evaluating 

n t iKrlV*) 

J[p(r|V0U),  p'rlVp)]  = log  - — Q- 

° ° U(r|V0°) 

where  the  observer's  state  of  knowledge  about  the  garget's  position  r is 
given  by  the  corresponding  rule  of  Bayes, 

t)  - ul?'vo0)  . 

0 J<£  U(p|*„°)  P{»gt|?) 

To  carry  out  such  a program,  one  must  possess  a reasonable  model  of  the  acoustic 
process,  capable  to  provide  the  range  and  bearing  dependence  of  the  conditional 
probability  density  p(Vq  |r)  . We  are  not  prepared  at  this  time  to  pursue 

this  problem  beyond  the  feasibility  remarks  above,  but  will  return  to  it  upon 
completion  of  a few,  somewhat  sampler,  investigations  that  naturally  suggest 
themselves  within  an  information-theoretic  approach  to  passive  acoustic 
detection.  It  is  important  to  notice,  however,  that  if  one  is  willing  to  view 
passive  acoustic  search  as  an  information-gathering  process,  there  exists  an 

I I 

entirely  natural  ard  conceptually  straightforward  way  of. positing  the  problem 
of  ranging  to  the  target.  In  fact,  estimating  the  target  position  is  in  no 
way  different  from  estimating  the  target's  existence  parameters.  As  we  have 
already  pointed  out  before,  this  highly  desirable  feature  of  the  information- 
theoretic  approach  to  detection  is  entirely  lacking  in  the  conventional  theory.  I 
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Indeed,  modern  passive  detection  apparatus  is  conspicuously  incapable  of  provid- 
ing ranging  information  or,  a previously  detected  target,  although  that  informa- 
tion was  definitely  resent  in  the  voltage  history  used  in  evaluating  the 
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SECTION  IV 

THE  CONVENTIONAL  APPROACH  TO  DETECTION 


The  preceeditKi  development  has  been  fundamentally  predicated  upon  the  sug- 
gestion that  searching  for  a submarine  is  primarily  a problem  in  information- 
gathering. Correspondingly , we  have  found  it  quite  natural  to  represent  the 
evoluation  of  the  search  process  by  the  time-dependence  of  the  observer's  state 
of  knowledge  concerning  the  complete  set  of  tactically  relevant  parameters.  In 
general,  the  attendant  body  of  theory  appeared  not  only  to  allow  the  evaluation 
of  the  searcher's  state  of  knowledge  along  some  voltage  sample  path,  but  also 
suggested  the  optimal  electronic  procedure  to  be  used  in  gathering  the  desired 
information. 

Notwithstanding  the  general  conceptual  appeal  of  this  approach  to  the 
acoustic  search  problem,  history  has  apparently  chosen  to  ignore  it.  To  be  sure  . 
many  of  the  concepts  presented  here  have,  in  some  form  or  another,  been  avail- 
able within  the  boundaries  of  conventional  theory;  but,  to  achieve  from  this 
vantage  point  comnlote  contact  with  the  classical  approach,  we  must  make  drastic 
and  often  quite  untenable  assumptions  the  impact  of  which,  while  hard  to  evaluate 
in  the  absence  of  tin  correct  result,  is  nevertheless  bound  to  be  significant. 

Within  the  i.or-pntiona  I theory,  a good  deal  of  relevant,  information  actually 
present  in  the  voltage  history  is  disregarded  in  that  the  complete  set  of  tactical 
parameters  is  generally  reduced  to  the  relatively  uninteresting  existence  parameter. 
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Despite  the  fact  that  the  observer  is  ultimately  interested  in  learning  the 
position  of  his  ta*qer,  it  appears  that  the  only  question  the  observer  is  in- 
structed to  ask  ot  tie  voltage  history  is  whether  or  not  a target  is  actually 
preent  n the  sut ruunaing  ocean.  To  describe  his  knowledge  concerning  the 
existence  parameter,  the  conventional  theory  insists  (reference  16)  that  the 
observer  make  objective  use  of  the  voltage  data  he  has  collected.  It  is  there- 
fore not  acceptable  to  refer  his  current  state  of  knowledge  to  an  a priori  one, 
unless  such  is  expressible  as  an  objective  frequency  statement  about  the  param- 
eter values  at  some  initial  time.  Correspondingly,  the  conventional  analyst, 
instead  of  using  nis  measurement  to  evaluate  the  amount  of  information  concern- 
ing the  presence  of  the  submarine,  would  use  his  measurement  to  evaluate  a 
related,  but  not  equivalent,  quantity  that  makes  no  reference  whatsoever  to  the 
observer's  a priori  state  of  knowledge.  This  quantity,  generally  referred  to 
as  the  likelihood  ratio  of  the  measured  sample,  no  longer  admits  of  any  simple 
interpretation  in  terms  of  concepts  relevant  to  the  original  search  question 
and,  therefore,  the  measured  value  of  the  likelihood  ratio  does  not  naturally 
relate  to  any  description  the  observer  might  consider  reasonable  for  his  know- 

I 

ledge  concerning  the  existence  parameter. 

lo  learn  about,  the  presence  of  a target  from  his  measurement,  the  conven- 
tional observer  must  a bitraril.y  correspond  the  value  of  the  likelihood  ratio 
to  the  possible  vihjes  of  the  existence  parameter.  In  general,  he  assumes  that 
the  target  is  presen;  if,  and  only  if,  the  likelihood  ratio  has  exceeded  some 
pre -ostatliS'he  I va  i o , the  magnitude  of  which  is  usually  so  determined  that  noise 


alone  wool  I seldom  cause  the  measured  value  to  go  above  it  (reference  17).  Th  i 
stochastic,  thresi  > id-ct ass iny  event  will  be  called  detection.  Needless  to  say 


the  conventional  way  establish  the  presence  of  a submarine  in  the  ocean 
implies  the  destructi  i of  some  finite  amount  of  information.  It  ir  therefore 
to  be  expected  that  a processor  operating  in  accordance  with  the  conventional 
framework  will  not  perform  the  task  of  determining  the  pre  ence  of  targets  as 
well  as  it  might  have  without  the  limitation  imposed  upon  it  by  the  information 
thresholding  procos. 

To  quantify  this  performance,  tne  conventional  theory  evaluates  the  proba- 
bility that  a detection  event  will  have  occurred  if  a target  were  actually  presert, 
and  calls  this  measu  e of  effectiveness  the  probability  of  detection.  We  ray 
illustrate  these  concepts  by  returning  to  the  simple  example  developed  in  our 
discussion  of  a deterministic  signal  where  a known,  undistorted  signal  is  mixed 
into  stationary,  white  Gaussian  noise.  Under  these  circumstances,  the  likelihood 
ratio  of  the  sample  v , 


whore  t represent  ..he  time-lengih  of  sample  viewed,  and  a detection  evont  ir 
caileo  if  q(t)  c^ce-ds  a predetermines  threshold  q-j  • 


In  principle,  there  are  a variety  of  ways  in  which  the  threshold  might 


meaningfully  be  chosen,  but  the  most  popular  practice  is  to  establish  a 
that,  will  maximize  tl  detection  probability  subject  to  a specified,  operationally 
acceptable  value  fo,  tliv  probability  of  false  alarm.  It  is  not  difficult  to 
determine  the  probability  density  function  of  the  correlation  integral  and, 
hence,  the  detection  probability  within  this  model.  In  fact,  since  the  measured 
voltage  is  assumed  to  be  a Gaussian  process,  the  correlation  integral  for  any 
fixed  value  cf  t is  a Gaussian  random  variable  of  mean  value 

(q(t))  ^ / dr  /v(t)\  s(t)  = f di  s2(t)  , 

\ /YES  O'  \ /yes  (T 

if  the  target  is  inueod  present,  and  of  corresponding  variance 


C 


) S ( T-j  ) S(t2) 


= j Nq  f dr  s2(t) 
c 0 


Therefore, 


1 \ -(q(t)-E)  /nqe 

p(qCt),  Yi  S';  ■*  ( VN0E  ) e 
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and  the  probability  oi  detection. 


Q 


d 


I dq  p(q  S YES) 


will  be  given  by 


where  E represents  the  energy  dissipated  during  the  observation  interval  in  a 
resistor  of  unit  resistance  when  a voltage  s(t)  is  put  across  it,  and  the 
threshold  value  qf)  is  related  to  the  false  alarm  probabil  ity, 

Qfa  2 I dq  p(q | NO)  , 

qo 


by  the  eqr  t'Cn 


•FA 


-erf 


Introducing  the  signal  to-noise  ratio 


p 


(t) 


,’F/U 


0 
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we  may  write  t.he  Jete:. 


.1  probability  as 


od(-) 


D-P;  t ) 


05) 


where  D is  related  lu  the  preassigned  false  alarm  probability  through 


1-erf  ( • ) 

\/2  I 


(16) 


It  is  to  be  noted  that  from  equation  (18)  the  probability  of  detection  for  a 
given  false  alarm  probability  is  a function  of  the  signal-to-noise  ratio  alone. 
This  property  continues  to  hold  true  for  an  entire  class  of  models,  but  such 
need  not  be  the  case  in  general. 


We  shall  close  Lb ' , section  of  our  review  with  yet  another  example  of  how 
the  conventional  thoc> / ascertains  the  presence  of  acoustic  targets.  Although 
this  model  has  little,  if  any,  intrinsic  claim  to  attention,  it  has  received 
sufficiently  widespreaj  application  to  warrant  inclusion  here.  Jt  assumes 
(reference  11)  that  the  ignal  is  represented  by  a white,  Gaussian  process  of 
mean  power  density  , and  that,  when  present,  this  signal  is  imbedded  in 
v.hi^e,  Gaussian  nm  ^ power  density  Nq  . Therefore, 
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2 Wt 

, ,ut  -,/2N0  £ “/ 

P(J|NO)=(?-U  e 


2Wt  2 

-|/2(vso>  S vj 


P{ ; l v es  ) = ( 2- m-fsry ) e 


and  the  likelihood  ratio  is  a monotonic  function  of  the  statistic 


2Wt  f t 

G = 1/Nfl  £ v 2 = 1/2N  dx  v2(t)  - 

D j=i  J 0J 


Under  the  null  hypothesis,  G is  x2-distributed 


p(G|N0)  * -TJFJ 

2 r(Ht) 


-G/2  nWt-l 
e b 


and,  correspondingly 


dG  p ( G | NO) 


» ^ f.  . ■.**■**•  ,k - * * ■■ 


Under  the  alternate  hypothesis,  G is  no  longer  x2  distributed,  but  the 
related  statistic 


G 


N, 


0 

+ s 


u 


G 


still  is.  Hence, 


p(G | YES)  - 


1 

Wt 

2 r(wt) 


e-G/2  GWt-l  dG 
dG 


and, 


Qd  = f dG  p(G  |YES)  . 

G0 

If  a variate  is  x2  distributed  with  degrees  of  freedom  n larger  than 
30,  the  square  root  of  that  variate  is  approximately  normal  with  mean  /2n  - 1 
and  unit  variance.  Therefore,  as  soon  as  2WT»30  , 


U 


FA 


f“d€  e-(P.-/4WFT)2/2 

/Mo 


C1.  = — 

d nr 

v 2ti  ~x. 


df;  e 


-U-/4WF1)2/2 
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which  corresponds  to  equations  (15)  and  (16)  if  we  replace 


p2(t)  4> 


(17) 


Naturally,  the  optimum  receiver  against  such  a signal  is  an  energy  integrator 
representing  at  the  output  terminals  the  statistic 


G(l)  - / oi  v:'(t) 

0 
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THE  SIGNA_  EXCESS  MODEL  OF  DETECTION 

The  examples  1"  ‘ofly  studied  above  appear  quite  useful  for  illustrating 
the  concept  . involved  in  the  conventional  theory,  but  can  hardly  be  considered 
representative  01  Ear  real  detection  process  for,  as  previously  indicated,  the 
signal  is  fundamental stochastic,  with  properties  that  are  strongly  related 
to  tne  physical  pm  r jmena  occurring  in  the  ocean.  A natural  attempt  at  accom- 
modating this  randomness  would  be  to  reformulate  in  the  conventional  language 
the  narrowuand  stochastic  signal  model  with  white  Gaussian  noise,  or  whatever 
other  model  migiu  be  suggested  by  an  analysis  of  ocean-fluctuation  phenomena 
(reference  lb). 

Practitioners  of  the  conventional  approach  have  found  it,  however,  more 
expedient  Lo  address  the  question  of  randomness  in  the  signal  from  an  entirely 
different  point,  c view.  According  to  this  "signal  excess  model"  of  acoustic 
detection,  it  is  the  randomness  in  the  signal-to-noise  ratio  process  rather 
than  the  randomness  in  the  signal  and  noise  processes  that  must  be  related  to 
the  parameters  characterizing  the  ocean  environment.  The  conventional  philosophy 
would  thus  instruc;  as  to  evaluate  the  mean  likelihood  ratio 

) - j f(  c)  A(  p)  > 

where  P(  p)  represents  the  probability  density  function  for  the  signal-to-noise 
ratio,  and  subsequently  to  compare  it  to  a preassigned  threshold  value  AQ  . 

The  probability  far  a detection  event,  to  occur  can  therefore  be  related  to  the 
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* 


measure  of  effectiveness  for  a processor  designed  to  perform  against  a given, 
non-random,  signal -to- noise  ratio. 


= J do  P{p)  Pr(A(p)  >AQ) 


If  is  important  to  nr t. ice,  however,  that  the  probability  for  a(p)  - Aq 
not  strictly  equal  to  ( p)  as  evaluated  in  the  previous  section,  because  the 
threshold  setting  is  now  chosen  to  maximize  P^  at  an  acceptable  level  of  false 
alarm  probability  rather  than  Qd(p)  itself.  With  that  understanding  in  mind, 
we  wri te 


P,  = f dp  P(p)  Q.(p)  , 

O'  u 

and  the  observer  ha^  agreed  to  be  satisfied  with  a threshold-setting  strategy 
that  performs  well  against  an  ensemble  of  possible  values  of  the  parameter 

distributed  according  to  P(p) 

Rather  than  use  any  of  the  models  for  Qrf(p)  discussed  previously,  it  is 

tht  oust on  to  take 


k 1 P > P !5 

Q .( n)  . 

i otherwise 


) 
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where  , , iv  tint  , 1 . 1- 1 -noise  ratio  at  which  the  probability  of  detection 

Q(J(  pi  becomes  eoual  t.u  h . Consequently, 

w» 

p„  * f it  ‘v  ; ■ 

r 

for  t,  r o/i,  . It  i,  further  customary  to  introduce  the  signal  excess  variable 
SE  = IP  log  5 


and  write 


Equation  (lb)  is  a succinct  statement  of  the  signal  excess  detection  model. 

Tne  signal  excess  variable  is  generally  written  as 

hi  - iO  j •.  , CD 

Mi  •• » r * 


RP  10  loq  r, 

■ *>  4 '(i •:  differential  (reference  19).  If  we  were  to  choose 

e.'i,.  i;,  actice,  the  white  Gaussian  signal  model  to  describe 

. ft-.*  *c  f.unitioii  differential  would  be  given  by 
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R1J  = 5 log 


DW 


with  the  preassigned  false  alarm  probability  related  to  D via  equation  (if,). 
The  first  term  in  t»e  signal  excess  equation  can  be  evaluated  from  the  sonar 
equation. 


THE  SIGNAL  EXCESS  PROCESS 

The  fundamental  ingredient  to  the  signal  excess  model  of  equation  (18) 
is  the  probability  density  function,  P(SE)  , describing  the  statistical 
properties  of  the  signal  excess  process,  SE ( t)  . Although  this  function 
should  be  specified  as  soon  as  some  acceptable  model  of  the  ocean  acoustic 
fluctuations  is  made  ivailable,  the  complexity  of  the  problem  and  a variety 
of  other  circumstances  have  encouraged  the  practice  whereby,  guided  by  the 
principle  of  simplicity,  we  choose  a mathematical  model  for  P(SE)  with  a 
few  undetermined  nanometers  and  subsequently  proceed  to  match  the  available 
degrees  of  freedom  to  the  data. 

The  Uncorrelated  Gm  .si  n Model 

A simple  and  h< nee  commonly  employed  model  chosen  for  P(SE)  is  a Gaussian 
model  with  mean  vlu  (SE(t))  as  provided  by  the  sonar  equation 

(sC(t>)  <F0M)  - (PL(t)) 


and  with  .idjustafi. 


lance 


MSL)  e-fSE-  (sE(t)))  /2 


2 2 
o 


> ?nn 


In  this  model,  the  s,\,ua)  excess  distribution  at  time  t is  independent  of  the 
previously  measured  values  for  the  signal  excess  and,  hence,  the  random  process 
SE(t)  is  uncorrelated.  The  detection  probability  is  trivially  evaluated  from 
equation  (21)  to  be 


, . erf  /.  imiLy  . 

Notice  that  (sE(t)\  = 0 corresponds  to  a probability  of  detection  of  ~ ; 

this  accounts  for  the  widespread  use  of  the  range  at  which  the  propagation  loss 
cun/e  crosses  the  mean  figure-of-merit  line  as  a measure  of  effectiveness  for 
sonar  performance. 

Values  for  the  variance  can  be  obtained  by  fitting  the  theory  to  specified 
,.cts  of  exercise  data,  but  due  to  the  general  paucity  of  reliable  data,  the 
choice  of  a is  mostly  a matter  of  taste. 

from  a tactical  viewpoint,  the  instantaneous  probability  of  detection, 
o (t)  , is  of  limited  value.  A widely  accepted  measure  of  tactical  effective- 

ness is  the  cumulative  probability  of  detection  (reference  20)  through  some  time 
interval  (0,t) 
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F(  t ) P*  ,f  [sE(t)  i.  0 , for  some  re  (0,t)] 


= 1 - | roi)  [SL(-i)  ••  0 , for  dll  re(0,t)]  . 

In  the  uncorrelated  Gaussian  model  discussed  here,  the  cumulative  probability  of 
detection  can  be  evaluated  as 


F«-  ’ - n ('-Vh» 

where  the  time  interval  (0,t)  has  been  partitioned  into  M equal  pieces,  and 
P^tj)  is  the  instantaneous  probability  of  detection  at  the  beginning  of  the 
jth  interval.  Naturally,  for  PcJ( t ) a continuous  function  of  time, 


Mm  Fm  - 1 

pl  oo 

it  0 
M.'.t  - t 

T he_  A - u J , 1 1 1 iu_  Mo de  I 

The  signal  excess  process  is  again  assumed  to  have  a normally  distributed 

/ \ 2 

marginal  with  mean  vSE(t)>  and  variance  o , but  the  signal  excess  values  £ 
different  time  instances  are  no  longer  statistically  independent  (reference  2:) 
In  fact,  the  process  is  fully  correlated  for  intervals  of  time  of  random  durati 
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distributed  Poisson  with  some  freely  adjustable  intensity  \ , and  fully  unc.or- 

related  otherwise. 

Ic  is  not.  dif.'i  -t.lt  to  evaluate  the  correlation  function  for  this  model 
of  signal  excess  variability.  Indeed  (reference  21),  let  us  consider  the 
statistical  dependence  between  the  signal  excess  values  at  t and  t + t , 
as  measured  by  the  correlation  function 

b(t,t+f)  = (SE(t)  SE(t+t))  - (SE(t))  (SE(t+t))  . 

With  probability  c‘AT  , there  will  have  been  no  jump  during  the  time  interval 
1 and,  hence,  with  that  probability 

SE(t+t)  = SE(t)  • 

With  probability  (l-e"AT)  at  least  one  jump  has  occurred  during  the  time 
interval  i , thus  totally  decorrelating  the  signal  excess  process 

Prob(SEft)  , SE(t+i))  = Prob(SE(t))  Prob(sE(t+T))  . 


Therefore , 
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<^Sf.(t)  SE(i.+t)\  - e >T  I °°  dx  x Prob(x) 

— co  * 

/ v 00 

+ (l-e  Al)  I dx  dy  xy  Prob(x)  Prob(y) 

— co 

= e~XT  (o2  + (sE(t))2) 

+ (l-e"AT)(sE(t))  (sE(t+x))  , 

and. 


b(t,t+T)  = o'  e'XT  + (sE(t))  ((SE(t+T)>  - (sE(t)))  e~XT  . 

\ 

If  the  mean  signal  excess  is  time-independent,  the  ^ - o jump  model  correlation 
function  becomes 


b { t ) - o?  e Vt  . 

By  definition,  the  instantaneous  probability  of  detection  is  the  same  here 
as  it  was  in  the  uncorrelated  model;  the  cumulative  probability  of  detection, 
however,  is  different,  and  no  longer  trivial.  Thus,  if  the  mean  signal  excess 
is  constant  in  time, 
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F(t)  1 Prob  (sE(t)  < 0,  for  all  -re(0,t)j 

= 1 * X,  Prob  (sE(t)<0,  for  all  Te(0,t)|j  jumps)  Prob  (j  jumps’) 

; - n > 

•.)  J 


so  that, 


F ( t ) = 1 - (!-rd  ) e"APdt  . 


From  this  equation  it  is  seen  that  the  mean  time  to  detection, 


One  can  similarly  show  that  if  contact  is  held  at  time  t = 0 , i.e.,  SE ( 0 ) >0 
then  the  probability  distribution  of  the  time  0 at  which  contact  is  lost  for 
the  first  time  is  given  by 


Prob  (e  < t)  = 1 - Pd  e"x^1-Pd^t  . 


The  next  level  of  complexity  of  known  results  concerns  the  case  when  the 
mean  signal  excess  is  time-dependent  and  unimodal,  that  is  SE ( t ) is  non- 
decreasing on  an  interval  (0,t0)  and  non-increasing  on  ( tQ ,«»)  . For  this 

case,  the  cumulative  probability  of  detection  can  again  be  computed  exactly 
(references  21  and  22), 


f J dT  pd(T> 

I'  - e 

| 

F(t)  " J 

-t 

-X  / di  Pd(x) 

J ' ( l*Pd^t0^ ) e 


for  arbitrary,  non-unimodal , mean  signal  excess  histories,  there  exists  no  closed- 
form  solution  to  the  problem  of  evaluating  the  cumulative  probability  of  detection, 
but  we  can  approximate  it  to  within  any  desired  degree  of  accuracy  (reference  23). 

lh_e_  Ga u ss- Ma rkov  Mode  1_ 

In  this  mode i , the  signal  excess  process  is  assumed  to  be  a Markov  process 
for  which  all  finite  dimensional  distributions  are  normal.  Oust  like  the  x - o 
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jump  model,  the  Gauss -Markov  model  describes  a correlated  signal  excess  process 
with  correlation  function  (reference  24), 


b(-r) 


0 


= o 


However,  unlike  the  jump  process,  there  are  no  useful  results  to  be  obtained  in 
closed  form.  The  principal  means  for  computing  the  cumulative  probability  of 
detection  are  approximation  methods  extensively  covered  elsewhere  (reference  22; 
and,  therefore,  we  shall  not  pursue  this  point  any  further.  We  will  note,  however, 
that  application  of  this  model  to  some  specified  circumstances  seems  to  indicate 
that  the  cumulative  probability  of  detection  for  the  Gauss-Markov  process  is 
larger  than  the  value  obtained  when  the  X - a jump  model  with  the  same  parameters 
is  used  under  the  same  circumstances  ( reference  25  ) . Furthermore,  one  shows 
(reference  26)  t.ha*  the  Gauss-Markov  process  allows,  on  the  average,  much  shorter 
contact  times  than  the  A - a jump  process.  In  fact,  for  large  positive  values 
of  (s£)  , the  mean  time  to  loss  of  contact  within  a Gauss-Markov  model, 


"G 
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j=! 


_ (<SE)/o>3 

2i/2  jr(iii) 


becomes 
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where  a,  represents  the  mean  time  to  loss  of  contact  within  the  X - a jump 

U 

model . 


The  Locj  Rice  Model 

All  stochastic  processes  discussed  above  have  been  so  designed  as  to  have 
a normal  marginal  probability  density  function  P(SE)  since  a large  body  of 
uata  seems  to  indicate  that  such  is  indeed  the  case.  Recent  studies  (reference 
27)  of  the  multipath  arrival  mechanism  for  randomness  in  the  received  signal 
seem  to  indicate,  however,  that  under  the  appropriate  circumstances,  the  signal 
excess  process  is  well  represented  by  (reference  28) 


SE( t)  - log  | ( P+x( t) )2  + y2(t) 

where  P is  an  arbitrary  non-negative  constant  and  x(t)  , y(t)  are  stationaiy 
Gauss-Markov  processes.  Although  pot  much  work  has  been  done  on  this  model  as 
yet,  one  can  show  (reference  29)  that  the  correlation  function  for  the  process  i' 

given  by. 


represents  the  correlation  coefficient  for  the  underlying 


CONCLUSIONS 


We  have  argued  that  acoustic  apparatus  should  be  employed  for  collecting 
information  about  a complete  set  of  tactically  relevant  parameters,  rather  than 
as  decision-making  devices.  Within  such  a framework,  information  that  is  gen- 
erally disregarded  by  the  conventional  procedure  becomes  available  and  can  be 
actively  used  in  whatever  decision  process  one  wants  to  adjoin  to  the  information- 
gathering stage  of  the  search  game.  This  additional  amount  of  information 
requires,  however,  that  one  possess  a physical  description  of  underwater  acoustics 
detailed  enough  to  provide  the  probability  distribution  for  the  acoustic  dis- 
turbance process  in  terms  of  all  tactically  relevant  parameters.  The  conventional 
theory  is  far  less  demanding  on  this  point  for  it  only  requires  descriptions  of 
the  acoustic  fluctuation  process  in  terms  of  ad  hoc  parameters  such  as  o and 
A to  be  fitted  to  the  available  data.  In  fact,  it  is  precisely  because  these 
parameters  are  not  directly  related  to  tactically  relevant  descriptions,  that 
the  conventional  approach  is  not  capable  of  obtaining  the  entire  amount  of 
information  available  in  the  voltage  process. 

It  is  not  difficult,  in  principle  at  least,  to  show  that  the  information- 
theoretic  approach  to  acoustic  search  does  indeed  lead  to  advantages  over  the 
conventional  framework.  To  do  so,  consider,  for  instance,  that  the  voltage 
sample  path  \l^  , has  been  drawn  from  an  ensemble  of  paths  characterized  by 
e=  YES  . An  observer  that  subscribes  to  the  information-theoretic  approach 


to  passive  acoustic  detection  would  use  his  processor  to  map  this  voltage  hist.o 
into  a number  describiig  his  state  of  knowledge  concerning,  say  e = YES  , in 
accordance  with  equation  (1), 


p(YES|V0°)  p(V0t| YES) 

p ( YES | VQ°)  p(V0t|YES)  + p(N0|V0°)  p^lNO) 


As  the  sample  path  V grows  in  length,  the  observer's  average  state  of  know- 
ledge (p(YES|VQt)  approaches  certainty  and  the  variance  of  his  sample-path 
state  of  knowleoge  about  this  mean  approaches  zero.  Therefore,  by  Chebyshev's 
inequality, 


Prob 


JplYESiV*)  - <^P  ( YES  | Vq1"  )>  )<e 


> 1 - 


Var(p(YES!V0t)) 


the  probability  that  the  sample-path  state  of  knowledge  p( YES  JVq11  ) approach 
certainty  with  increasing  viewing  time  gets  arbitrarily  close  to  unity.  The 
convent ir  ;al  cbser  c.  • r,  the  other  hand,  will  employ  his  processor  to  map  the 
voltage  sample -pa  a > to  a succession  of  YESes  and  NOs  by  automatically 
thresholding  a runnv  -window  statistic.  Thereupon,  depending  upon  the  rela- 
tive frequency  of  Er:  in  the  succession,  he  will  subjectively  translate  this 

success  irv  :nt."  a cl  imate  decision  regarding  the  presence  of  the  target. 


To  make  the  corns*-  i son  between  the  2 observers  meaningful,  it  seems 
reasonable  to  assume  a at  the  conventional  oDserver  will  use  the  rule  of  Bayes 
ro  t cans la  the  sac  sion  of  YESes  and  NOs  he  has  available  into  a state 


of  knowledge  concerning  the  existence  parameter.  Thus, 


p ( YES  j Dq1)  = 


p(YES|D0°)  pfD^IYES) 

p(YES|D0°)  p(D0VES)  + P(N0|D0°)  p(D0t|N0) 


where  represents  a sample-path  succession  of  ones  and  zeroes  corresponding 

to  the  succession  of  YESes  and  NOs  generated  by  automatic  thresholding,  and, 
for  a running  window  of  size  W , 


P(D| YES)  = / dv  p(v|YES)  IT  0 

J=1 

is  the  probability  of  a discrete  realization  D of  the  sample-path  Dq^  . 

It  is  not  our  intention  to  accomplish  the  numerical  comparison  between 

these  two  processing  schemes  at  the  present  time,  but  it  should  be  obvious  by 

t 

now  that  the  only  difference  between  them  lies  in  the  replacement  of  p(Vq  ] f) 
with  p ( Dq ^ | t)  . As  equation  (19)  indicates,  the  conventional  processor 
performs  an  averaging  operation  over  all  those  voltage  sample-path  Vq*  that 
produce  the  same  sample-sequence  before  evaluating  the  state  of  knowledge 

p(YES(Drjt)  , and  therefore  neglects  a certain  amount  of  statistical  detail 
available  in  the  voltage  history  it  reads. 

If  we  accept  the  information-theoretic  approach  to  the  submarine  search 
problem,  the  conventional  concept  of  a detection  event,  together  with  most  of 
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the  analytic  modeling  developed  to  assess  the  performance  of  sonar  gear,  should 
be  suitably  generalized.  Consequently,  much  work  remains  to  be  done.  To  begin 
with,  a thorough  physical  description  of  ocean  phenomena  relevant  to  acoustic 
propagation  should  be  attempted.  Dr.  William  Hurley  has  recently  been  assigned 
to  participate  in  the  Acoustic  Detection  Project  and  plans  to  investigate  the 
feasibility  of  methods  ether  than  the  wave  equation  for  connecting  the  statistical 
properties  of  the  ocean  and  its  boundaries  to  the  statistical  properties  of  under- 
water acoustics.  Next,  the  information  process  itself,  given  a probability  dis- 
tribution for  the  voltage  history  received,  must  be  addressed  in  some  detail  to 
carry  through  the  numerical  comparison  indicated  above.  Finally,  we  envision 
utilizing  this  framework  in  conjunction  with  some  acceptable  decision-making 
scheme  to  deal  with  certain  of  the  Navy’s  major  tactical  necessities  such  as 
passive  ranging  to  a target. 
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